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ABSTRACT 



i-C Aims. Given that in most cases just thermal pressure is taken into account in the hydrostatic equilibrium 
equation to estimate galaxy cluster mass, the main purpose of this paper is to consider the contribution of all 
q three non-thermal components to total mass measurements. The non-thermal pressure is composed by cosmic 
5— l ■ rays, turbulence and magnetic pressures. 

t/3 Methods. To estimate the thermal pressure we used public XMM-TVewfon archival data of five Abell clusters to 
C3 , derive temperature and density profiles. To describe the magnetic pressure, we assume a radial distribution for 
1 the magnetic field, B(r) oc pg, to seek generality we assume a within the range of 0.5 to 0.9, as indicated by 
observations and numerical simulations. Turbulent motions and bulk velocities add a turbulent pressure, which is 
■ considered using an estimate from numerical simulations. For this component, we assume an isotropic pressure, 
("^ fturb = |p g (af + af). We also consider the contribution of cosmic ray pressure, P cr oc r -0,5 . Thus, besides the 
gas (thermal) pressure, we include these three non-thermal components in the magnetohydrostatic equilibrium 
\Q • equation and compare the total mass estimates with the values obtained without them. 

CD Results. A consistent description for the non-thermal component could yield a variation in mass estimates that 
. extends from 10% to ~30%. We verified that in the inner parts of cool core clusters the cosmic ray component 
i—( ' is comparable to the magnetic pressure, while in non-cool core clusters the cosmic ray component is dominant. 
^\ t For cool core clusters the magnetic pressure is the dominant component, contributing more than 50% of the 
CD ■ total mass variation due to non-thermal pressure components. However, for non-cool core clusters, the major 
influence comes from the cosmic ray pressure that accounts for more than 80% of the total mass variation due 
to non-thermal pressure effects. For our sample, the maximum influence of the turbulent component to the total 
mass variation can be almost 20%. Although all of the assumptions agree with previous works, it is important 
' to notice that our results rely on the specific parametrization adopted in this work. We show that this analysis 
can be regarded as a starting point for a more detailed and refined exploration of the influence of non-thermal 
pressure in the intra-cluster medium (ICM). 



1. Introduction 

Clusters of galaxies are powerful tools for investi- 
gations of cosmological interests. The evolution of 
the mass function of a cluster is highly sensitive to 
cosmological models since the matter density con - 
trols the rate at which structures grow (|Voitll2005( ). 
In order to use clusters of galaxies as observational 
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probes of dark energy in the Universe and to in- 
vestigate the structure formation history including 
baryonic hydrodynamics, the non-thermal contri- 
bution must be well understood and quantified. 

X-ray data are one of the most often methods 
used to determine the mass distribution of clus- 
ters of galaxies. To do so, hydrostatic equilibrium is 
usually assumed, and the observed gas density and 
temperature profiles are used to compute the ther- 
mal pressure. In most cases, only the gas (thermal) 
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pressure is considered to evalu ate the dynamical 



masses of galaxy 


clusters (e.g., David et al. 


1995; 


White fc Fabianl 


1995; Finoeuenov et all 


2001; 


Reiprich & Bohringer 


2002). However, there is also 



a non-thermal pressure (-Pnt), composed by the 
magnetic (Pb), turbulent (Pturb) and cosmic ray 
(Per) components that is frequently assumed to be 
negligible and thus ignored. As a consequence of 
this, today the accuracy of the hydrostatic mass es- 
timates is limited by the non-thermal pressure from 
these components. 

Despite the difficulty to reliably calculate the 
small-scale properties of the magnetic field, the ex- 
istence of intra-cluster magnetic fields is well estab- 
lished from the studies of the rotation measure of 
polarized radio frequencies and s ynchrotron emis- 
sion from diffuse sources (e.g .. | Andernach et al 

' [Tax 



1988; 
2002 : 
other 



Giovannini et al 



Govoni fc Feretti 




Taylor et alj 11994 



More recently, an- 
indication that the intra-cluster medium 
(ICM) is permeated by a magnetic field came from 
the studies of X-ray cold fronts (sharp discontinu- 
ities in X -ray surface brightness profile and tem- 
perature, iMarkevitch fc VikhliniiJl2007t ). In these 
cases, a parallel magnetic field can suppress trans- 
port processes in the ICM, making it difficult to 
mix different gas phases during a cluster merger. 
Even a very weak magnetic field can effectively in- 
hibit transport processes such as th ermal conduc- 
tion and the settling of heavy ions (jSarazi n 198G; 



Soker fc Saraziij |l990) 



Strong magnetic fields can make a signif- 
icant contribution to the gas pressure sup- 
contributing with a 
component in the magnetohydro- 



port IjLoeb fc Mao 
non-thermal 

static equilibr i um e quation (jDolag et alJ [2001b; 
IVogt fc Enfilml 120051 ). Indeed, magnetic fields as 
high as 10 - 100 txG were found in Hydra A 
dTavlor fc Perlevl Il993h . Cygnus A dDreher et al 



19871) an d m3C 295 (IPerlev fc Tavlorl Il99lh 
Moreover, Dolag et al. ( 19991 ) performed numerical 



simulations and found that even clusters with an 
overall small magnetic field can be penetrated par- 
tially by regions of high magnetic fields. 

Although on average the magnetic pressure 
in simulations is much smaller than the thermal 
pressure (~ 5%, iDolag fc Schindlerl l2000h . there 
are domains of high magnetic fields approaching 
or sometimes even exceeding equipar tition with 
the thermal energy. Previous studies (IDolag et al 



2001aHColafrancesco fc Giordanoll2007 . among oth- 
ers) have analyzed the effects of the magnetic pres- 
sure in simulated galaxy clusters. 



Magnetic fields and turbulence are possibly re- 
lated to one another. It seems plausible that the tur- 
bulent motions in the ICM can maintain the mag- 
netic field b y converting kinetic energy into mag- 
netic energy ( Sanchez- Salcedo et al.lll999h . The ob- 
served small-scale turbulence in the ICM can be 
due to bulk velocities and ongoing merger of sub- 
structures. Gas turbulence on small scales can also 
be driven directly by motions of galaxies, as for in- 
stance by jets and bubbles from the active galactic 
nuclei (AGN. IChurazov et al.l 120021 ) . although the 
latter may be confined to the inner regions of the 
cluster ( Lau et al. I l2009l ). The presence of random 
gas motion can also contribute to the pressure sup- 
port in clusters of galaxies. 

The chaotic nature of the ICM magnetic field 
makes it difficult for energetic particles to scape 
from the cluster, and thus cosmic-ray protons would 
be confined for timescales exceeding the Hubble 
time. The electron cosmic rays, on the other hand, 
have collisional and radiative lives much shorter. 
Thus, since the ICM is permeated by significant 
magnetic fields, one would expect the cosmic ray 
pressure to have some relevance in its support 
against gravity. 

To consider deviations from the standard as- 
sumptions in computing cluster total mass, the 
main aim of this work is to analyze the effects of 
non-thermal pressure, that is to take into account 
magnetic, turbulent and cosmic ray components. 
Hydrostatic masses were derived using X-ray obser- 
vational data for five Abell clusters: A496, A2050, 
A1689, A2667 and A2631. To do so, we use tem- 
peratu re and density pro file fits from a previous 
work ([Lagana et al . 2008) and wc introduce the 
Pnt contribution in the magnetohydrostatic equi- 
librium equation. For these five clusters, we com- 
pare masses determined considering non-thermal 
pressure (Mntp(^)) with their hydrostatic values 
(M(r)). 

The paper is organized as follows. We show the 
data sample in Sect. [51 The non-thermal compo- 
nents are described in Sect. [3] In this section we 
describe the structure of the intra-cluster magnetic 
field, the turbulence in the ICM and the cosmic ray 
component. In Sect. (U we present the method of 
determining the cluster mass, including the effects 
of the Pnt- Our results, as well as a discussion of 
them are presented in Sect. [5] and our conclusions 
in Sect. [6] 
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Table 1. General cluster properties. 



Cluster 


R.A 


DEC 


z 


rsoo 




(J2000) 


(J2000) 






A496 


04 33 37.1 


-13 14 46 


0.033 


1480 


A2050 


15 16 21.6 


+00 05 59 


0.1183 


2172 


A1689 


13 11 34.2 


-01 21 56 


0.1823 


1785 


A2667 


23 51 47.1 


-26 00 18 


0.23 


2153 


A2631 


23 37 39.7 


+00 17 37 


0.273 


1976 



2. Data sample 

The objects in our sample are within the redshift 
range of 0.03 < z < 0.3 and are drawn from a set 
of Abell clusters with available data in the XMM- 
Newton public arc hive. These clu s ters w ere previ- 
ously analyzed by lLagana et alj (|2008l ). who de- 
rived the density profiles fit parameters to compute 
the total mass. Although we have not used the ob- 
ject morphology as a criterion for the cluster selec- 
tion, all these clusters except A2631 have appar- 
ently symmetric X-ray isophotes, suggesting that 
they are relatively relaxed. The deviations in the 
surface brightness profile of A2631, although clearly 
present, are not very large and do not invalidate the 
assumption of spherical symmetry. However, we are 
aware of the fact that it may affect total mass re- 
construction, ac counting for underestimated mass 
determinations (jPiffaretti fe Valdarninill2008l ). 

In Table [TJ we present the five Abell clusters 
used in this work, specifying rsoo, the radius inside 
which the mean density exceeds the critical den- 
sity by a factor of 500. All masses are computed 
inside rsoo, as it is the largest radius for which the 
current X-ray data re quire no model extrapolation 
(Vikhlinin et a . 2006) and is about the virial radius 



Lacev fc ColeT Tl993) 



Usually, the mass of a cluster is determined un- 
der the assumption of hydrostatic equilibrium with- 
out the contribution of non-thermal pressure. In 
this case, the total mass relies on the temperature 
and density profiles. 

Satellites with better spatial resolution 
(like the XMM-Newton and Chandra) showed 
a significant difference between the surface 
brightness profile data points and the /3-modcl 
( Cavaliere fe Fusco-Femianc] 1976 . 19781) a t small 



radii for cool core clusters 
Xue fc WiJ l2000h 
difference 



the 



Jones fc Formanlll984 : 
Based on this observational 
/3-model was used to describe 



the density distribution of non -cool core (NCC) 
clusters, while the S ersic model (|Pislar et al.lll997 : 
Demarco et al.ll2003r ) was used to characterize cool 



core (CC) clusters. For A2050 and A2631, the gas 
density (p g ) is described by 



Pg(r) = Po ( 1 + ^ 



-3/3/2 



(1) 



where po and r c are the central gas density and 
the gas core radius, respectively. The (3 parameter 
determines the power-law behavior at large radii. 
For A496, A1689 and A2667 (CC clusters) the gas 
density profiles were fitted by the Sersic model given 

by 



Pg( r ) 
where p' 



r 
1 a 



exp 



(2) 



p/2, y = l-0.6 097i/+0.05563i/ 2 and . 



a' 2 X I V (jDurret et al.ll2005t l. The best-fit parameters 
were determined from the X-ray surface bright ness 
profiles and were given by lLagana et al. ( 2008 ) . 



3. Non-thermal components 

In this section we describe each non-thermal com- 
ponent considered to contribute to the pressure sup- 
port. 



3.1. The magnetic profile 

In the 80s, iJaffe (|l980h suggested that the intra- 



cluster magnetic field distribution should depend 
on the thermal gas density and on the distribution 
of massive galaxies, which means it would decline 
with cluster radius. Cluster observations provided 
constraints on the radial gradient of the cluster 



magne t ic field (Brunc tti et al.M2001t iGovoni et al 



2001at iFeretti et all 120041 ). The intensity of the 
magnetic field was found to decrease smoothly with 
the cluster-centric radius, with a trend similar to 
that of the thermal gas. 

From magnetohydrodynamic (MHD) cosmolog- 
ical simulations, an important characterization of 
the cluster magnetic d istribution was made by 
Dolag et al] (|l999l . l2002[) . They studied the correla- 
tion between X-ray surface brightness and Faraday 
rotation measurements (FRMs) in clusters provided 
by X-ray and radio observations as well as from 
models for the ICM. These authors performed cos- 
mological MHD simulations in order to recover the 
correlation between these quantities. They found a 
relation between magnetic fields and the gas den- 
sity of the cluster, suggesting that the cluster mag- 
netic fields may span a wide range of spatial scales 
with a strength t hat decreases with d istance from 
the cluster center. iMurgia et al.l (|2004l) used numer- 
ical simulations to investigate the relation between 
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magnetic fields and Faraday rotation effects in clus- 
ters. These latter authors compared their simula- 
tions with polarization properties of extended clus- 
ter radio sources in radio galaxies and halos. They 
considered that the intensity of magnetic fields de- 
creases from the cluster center in agree ment with 
pr evious results (IDolag et al.lll99 9. 2002). 

Dolag et all (|l999h found that the observed 



intra-cluster magnetic field can be reproduced by 
the evolution of an initial magnetic field at red- 
shift f 5 that was amplified by compression during 
the cluster collapse. One of their important results 
was that the intra-cluster magnetic field strength 
is proportional to the gas density at any point 
(B (r) oc p(r)). 

IColafrancesco fe Giordano! (|2007l ) studied the 
influence of magnetic fields on the main structural 
properties of virialized groups and clusters, assum- 
ing that it scales with a density of B(r) oc p a (r), 
as previously proposed. The same power law de- 
pen dence on the density was used bv I Zhang ( 2004h 
and iKoch et all (|2003h to estimate the effect of 
the intra-cluster magnetic field on the Sunyaev- 
Zel'dovich power spectrum. 

Motivated by these previous works, we assumed 
a parametric form for the radial distribution of the 
magnetic field 



B{r) = B 



PgiT) 
Po 



(3) 



where Bo is the central value of the magnetic field 
and a is the shape parameter. Unfortunately, there 
are no measurements for the magnetic profile for 
any of the clusters in our sample. Thus, we have to 
use results from the literature to constrain B and 
a. 

The effective strength and structure of these 
magnetic fields provide the main challenge, because 
different methods of analysis give different values 
for magnetic strength. An estimated consideration 
of the equipartition of the magnetic field strength 
averaged over the entire halo v olume gives magnetic 
field strengths of ~ 0.1 - ljuG (jGovoni et alj|2001b : 
Murgia et al.l 12004, and references therein) 



Feretti et al. ( 1999al ) estimated that the mag- 



netic field in the ICM of A119 s hould range between 
5 - 10/iG. iBaechi et all (|l998l) found B w 1/iG for 
the cluster-sc ale magnetic fi eld st rength. In a more 



recent work, Clarke et al.l (|200ll ) studied a sam- 
ple of 16 "normal" low-redshift (z < 0.1) galaxy 
clusters, finding that the ICM is permeated with 
magne tic fields at levels of 4-8 [iG. iTavlor fc Perlevl 
(1993) found higher central val ues, B ~ 6 — 30^G, 
for the ICM magnetic fields. lAllen et all (|200lh 



claimed that the central value of the magnetic fields 
can be B — 12/iG. FMRs of radio sources pro- 
vide magnetic fields of ~ 5 — 30 jiG in cooling 
flow clusters (e.g., 3 C 295 dAllen et alfeOOll). Coma 
dFeretti et alJll995h and Hydra A (jTavlor fe Perlev 
1993)) where extremely h igh F RMs have been re 



vealed. ICarilli fc Tavlorl (|2002t ) affirmed that its 
strength in the center of cooling-core clusters can 
reach levels of 10 — 40/iG. 

On the other hand, lower magnetic fields (~ 2 — 
8fxG) have also b een detected in clusters without 
cooling flows (e.g.jFeretti et al.lll999bl Taylor et al 



120011 lEilek fc Oweij|200aT 

The magnetic field strengths obtained from 
FRMs arguments are higher than the values derived 
either from the radio data or from inverse Compton 
X-ray emission. The values deduced from radio syn- 
chrotron emission and from inverse Compton refer 
to averages over large volumes. Instead, FRMs es- 
timates give a weighted average of the field and gas 
density along the line of sight and could be sensi- 
tive to the presence of filamentary structure in the 
cluster. They could therefore be higher than the 
average cluster value. However, as pointed out by 
Carilli fc Taviorl (|2002l) . all of these techniques are 



based on several assumptions. For example, the ob- 
served FRMs have been interpreted until now in 
terms of simple analytical models which consider 
single-scale magnetic fields, while equipartition cal- 
culations in radio halos assume spatially uniform 
magnetic fields. 

There are not many works that studied the 
power spectrum of the intra-cluster m a gnetic field 
fluctuations. Howeve r Enfilin fc Vogt ( 20031 ) and 
IVogt fc Enfihnl (|2003l |2005|) by using a new semi- 
analytical technique showed that, for those cluster 
sources for which a very detailed FRM image is 
available, the magnetic field power spectrum can 
be estimated. 

To take into account all these observational re- 
sults, we considered conservative values in our anal- 
ysis, adopting for cooling core clusters values for B 
ranging from a low value of 5/iG up to a strong value 
of 30^G. But for non-cooling core clusters we con- 
sidered the central strength of the magnetic field to 
vary from 2/xG < B < 8^iG. The values adopted 
agree with both observations and with numerical 
simulations that considered higher values for mag- 
netic fields in the core of the cluster and a decrease 
towards the outskirts. 

If the ratio of magnetic and thermal pressure is 
constant throughout the cluster, then a = 0.5, if 
we have an homogeneous seed magnetic field com- 
pressed during cluster formation, one expects a — 
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2/3. IVoet fc Enfjlm] (|2005l ) performed a Bayesian 
maximum likelihood analysis of Faraday rotation 
measure in order to derive a power spectrum of clus- 
ter magnetic fields. In their analysis they used three 
different values for a = 0.1, 0.5, 1.0. They concluded 
that values of a = 1.0 are unlikely, but models 
with a — 0.1 — 0.5 match very well wi th their cal- 



culations. From observational results, Dolag et al 
(1200111 found a = 0.9 for A119 and a = 0.5 for 3C 
129. If one combines the measurements of the four 
clusters the data analysis performed by these latter 
authors leads to a slightly lower slope, of this cor- 
relation (a = 0.8). Therefore, to be conservative, 
we decide to vary the shape parameter within the 
interval of 0.5 < a < 0.9 to take into account all 
the above results. 

In Fig. Q] we show the mass profile variation due 
to the magnetic pressure as the only non-thermal 
component compared to the hydrostatic mass pro- 
file. From this Figure we see that the magneto- 
hydrostatic profiles for NCC clusters (A2050 and 
A12631) present little difference when compared to 
the hydrostatic profile. This little influence on the 
magnetic pressure in NCC cluster is due to the cen- 
tral value assumed for these clusters. As the central 
strength of the magnetic field is lower in NCC clus- 
ters, the influence of the magnetic pressure is also 
lower compared to a CC cluster. For CC clusters 
we note that the difference between the hydrostatic 
and the magnetohydrostatic profiles becomes more 
pr onounced for large radi i, that is for r > 0.5 rsoo- 
(|2008h 



Churazov et al 



measured the contribu- 
tion of the non-thermal pressure in two early-type 
galaxies that reside in the center of two nearby 
cool core clusters (Virgo and Fornax). With a sim- 
ilar approach, these authors considered the con- 
tribution of all non-thermal components in the 
same system, combining the contribution of cos- 
mic rays, magnetic fields and microturbulence to 
the total pressure. They suggested that the Pmt 
component can account for 10% of the gas thermal 
pressure in the core of these galaxies (NCC 1399 
and M87) . Assuming that the mag netic component 



is the only non-thermal pressure, IChurazov et al 
(2008) constrained the upper limits on magnetic 
field to be ~ 20-30 /xG. Besides. lAjeilo et all (|2009l ) 
reported the detection of ten merging-clusters in 
the 15-55 keV energy band using Swift's Burst 
Alert Telescope (BATJ). These authors coupled ra- 
dio synchrotron emission (inverse Compton) with 
X-ray data, putting the constraints on the lower 
limit of the magnetic field to be typically between 
0.1-0.5/iG. We can see from these Papers mentioned 




0.1 0.3 


0.5 

r / r 500 


0.7 


0.9 


A 1689 
















— ^Th 








0. 1 0.3 


0.5 


0.7 


0.9 




|http: / /swift .gsfc.nasa.gov/dQcs/swift/about_swift/bat_desc.htmll 



'I'm 

Fig. 1. Variation in the mass profile due to the mag- 
netic pressure as the only non-thermal component for 
the clusters in our sample. The red lines represent the 
hydrostatic mass profiles, while the green lines show the 
maximum variation of mass profiles due to the contribu- 
tion of the magnetic pressure. The blue zone represents 
the region of mass profile variation depending on the 
central value strength (-Bo) and the shape parameter a. 



above that both the upper and lower limits adopted 
in this work are in accordance with the recent lit- 
erature. 

3.2. Turbulence in clusters of galaxies 

With the advent of high-resolution observations 
from the Chandra and XMM-Newton satellites, 
temperature maps of the X-ray emitting gas 
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have shown that even apparently relaxed clusters 
could have small substruc tures as imprints of re- 
cent minor-mergers (e.g. 



Finoguenov et al.l 2005 



Durret fc Lima Netolf2008h 

It is widely accepted that the ICM is proba- 
bly turbulent, and mergers of galactic sub-clusters 
may be on e of the majo r energy injection mecha- 
nism s (see Sarazin 2002 : Brunetti 20031 iLazarianl 
20061 and references therein). Random gas mo- 
tions can also maintain and ampli fy cluster mag- 
netic f ields via dynamo processe s ( Roettiger et al.l 



In order to quantify the importance of pressure 
support from random gas motions in clusters, we 
can write the following relation for the isotropic 
turbulent pressure Pturb (|Lau et al. 2009t) : 



-Pturb = gPgK 2 +CT t 2 ), 



(5) 



19991: 

to the 



Subramanian et al. 2006) and contribute 
acceleration of co smic rays in the ICM 



(jBrunetti fc Lazar ian 2007). 

Generally the models assume a scenario with the 
scale for the injection of energy of 100-500 kpc and 
the injection velocity of the order of 10 3 km/s. 

Since the rate of dissipation of the turbulent en- 
ergy cannot exceed the X-ray luminosity of the clus- 
ter (Lx) in a steady-state, i.e., ^Vq/Iq ^ Lx/M g , 
where vq and Iq are the turbulent speed and the 
scale respectively, and M g is the gas mass, we have 
an upper limit on the tur bulent velocity as follows 



where ay and at are the radial and tangential 
dispersion velocity of the intra-cluster gas respec- 
tively. For 16 simulated clusters with virial masses 
within th e rang e of (5 x 10 13 - 2 x 10 15 ) M Q , 
lLau et ah ( 20091 ) found that gas motions contribute 
up to ~ 5% — 15% of the total pressure sup- 
port in relaxed clusters. Thus, on average the to- 
tal mass estimate is biased low by about 8 ± 2% 
(at rsoo) in relaxed systems and 11 ± 6% in un- 
relaxed systems. These results agree with previous 
studies (lEvrardll990HRasia et alj2004HNagai et al' 



( Subramanian et al. 20061 ) 



v < 180 



km / l 



L 



x 



(- 



s V 2 00 k Pc7 \W 45 erg/s) \10 1A M e J 'haloes with 



l2007t IPiffaretti fc Valdarninill2008l) . with contribu- 
tions increasing along the radius. 

Note that random gas pressure and its gradient 
is sensitive to small-scale clumps and any pressure 
inhomogeneity, and these sources could potentially 
bias the measurements of the pressure gradient and 
hence the hydrost atic mass estima te. In order to 
(minimize this bias Lau et al. ( 20091 ) removed sub- 



Norman fc Brvanl (|l999h found that the ICM 



becomes turbulent during cluster formation, with 
turbulent velocities of about 400 km/s within 1 Mpc 
from the center of a cluster and eddy sizes ranging 
fr om 50 to 500 kpc . In th e cluster merger model 
of Ricker fc Sarazin ( 200 if ), they found large-scale 
turbulence with eddy sizes up to several hundred 
kiloparsecs and turbulent velocities of ~ 100 — 
400 km/s. 

From the analysis of pressure fluctuations as 
revealed in X-ray observations, ISchuecker et al 



a mass greater than 10 h Mq and 
the mass within their tidal radius from their calcu- 
lation (see their paper for further details). 

In our analysis, we used radial and tangential 
velocity dispersion prof iles (see Fig|2l) ba sed on the 
numerical simulation of lLau et al. ( 20091 ). We show 
these profiles in Fig. [2] where the velocity 



.-.(..in 



GM 

^500 



(6) 



( 20041 ) argue that the integral turbulent scale in the 
Coma cluster is close to 100 kpc, and they assume 
a turbulent speed of 250 km/s at that scale. 

Although a number of other studies have exam- 
ined random gas m otions and their e ffect on the 



mass estimate fe.g.. IRasia et al . 2004; iDolag et al.. 
20051 : Rasia et al. 20061 ). most of them have used 



simulations with SPH gas dynamics. The magni- 
tude and effects of gas motions in such simulations 
depen ds on the specific treatment of artificial vis- 
cosit y (IDolag et al.ll2005l) . Thus, we used the results 
from lLau et alJ (|2009l ). which employed simulations 
with Eulerian gas dynamics with very low numer- 
ical viscosity. This approach is therefore useful in 
evaluating possible differences between numerical 
techniques and systematic theoretical uncertainties. 



is the circular velocity at rsoo- 

In Fig. [3] we show the mass profile variation due 
to the turbulent pressure as the only non-thermal 
component, compared to the hydrostatic mass pro- 
file. From this Figure we see that the influence of the 
turbulent pressure in the mass estimates is small re- 
gardless of weather the cluster is a non-cool core or 
cool core cluster. Our sample has masses within the 
range of (0.3—12) x 10 14 M Q , and we found that our 
mass estimates can be biased low by about ~5%. 

3.3. Cosmic ray pressure 

Cosmic ray protons can play an important role 
within ICM contributing to the equilibrium on the 
pressure support. Cosmic ray protons (CRp) and 
cosmic ray electrons (CRe) can be injected into 
the ICM by three different processes which produce 
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Fig. 2. Velocity d ispersion p r ofiles derived from numer- 
ical simulation of lLau et al.l ((2009) . Upper panel: tan- 
gential velocity profile used to estimate the turbulent 
pressure. Lower panel: radial velocity profile used to es- 
timate the turbulent pressure. 



morphologically different radio signature (jBrunetti 
20031: IPfrommer fe En6rin1l2004 : 



shock acceleration: natural acceleration 
mechanism such as structure formation 
and merger shocks dPtuskin et all 120091 : 



Ostrowski fc Siemieniec- Oziebloll2002ri 



reaccelerated electrons: reaccelerated rel- 
ativistic particles injected by sources like 
radio galaxies, supernova remnants, merger 



shocks, galactic w i nds, e tc. JEnsslin et al.lll997 



Berezinskv et al 
2004 . 



1997: Pfrommer fc Enfilin 



— particles of hadronic origin: CRp can eventually 
interact with the thermal ambient gas producing 
secondary electrons, neutrinos and 7-rays in an 
inelastic collision which generates a radio halo 
through synchrotron emission. 

Evidence of non-thermal electrons in the ICM 
exist in the form of synchrotron radio emissio n 
( Feretti fc Giovanninil l200i iBrunetti et "all l2008h . 
excess of ex treme-ultraviolet (EUV) and hard X-ray 
radiation dBowver fc Berghdferl 1 19981: iLieu et al 



19991: iFusco-Femianol 120041 ) Another consequence 
of the presence of cosmic rays in the ICM is the 
production of gamma rays. 

In order to consider the contribution of cosmic 
ray pr essure to the ICM, we followed the prescrip- 
tion of lAndo fc Nagai (|2008f) . We can define a rela- 
tive contribution of the cosmic ray pressure support 



Y„ 



p 




'I'm 

Fig. 3. Variation in the mass profile due to the tur- 
bulent pressure as the only non-thermal component for 
the clusters in our sample. The red line represents the 
hydrostatic mass profile, while the green line shows the 
maximum variation of the mass profile due to the con- 
tribution of the turbulent pressure. The blue zone repre- 
sents the region of the mass profile variation depending 
on the radial and tangential dispersion velocities. 



where this ratio can be parametrized using a power 
law 



Y P (r) =Y p0 [- 



(8) 



(7) 



where the s ubscript rep r esents values in the cen- 
tral region. Siiacki et al. |2008), using numerical 
simulations, followed the evolution of the cosmic- 
ray plasma inside the bubbles, considering both 
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its hydrodynamical interactions and the dissipation 
processes relevant to the cosmic ray population. 
They found that cosmic ray pressure is most rele- 
vant in the center of the clusters, being comparable 
to the thermal pressure for r < 50h~ 1 kpc. 

Recent studies have highlighted that CRs can 
be dynamically important in galaxy clusters be- 
cause they put constraints on the fraction of cosmic- 
ray pressure with respect to the thermal pressure 
(Y p ). Since the Energetic Gamma Ray Experiment 
Telescope (EGREtQ) did not detect 7-ray emissio n 
from clusters in the GeV band ( Reimer et alj|2003 ) , 
constraints on the fraction of cosmic ray pressure 
have been placed in cosmological simulations of the 
large scale structure. In nearby rich clusters, this 
component shoul d amount to about ~10-26% of 
thermal pressure fensslin et al. 1997 : Miniati et al 



200lUMmiatill2003t ). By comparing the integrated 
7-ray flux above 100 MeV t o EGRET upper limits, 
Pfrommer fe EnBlinl (|2004l ) constrained the CRp 
scaling parameter in their simulation of nearby 
cooling-flow clusters. Thus they were able to infer 
that the P CI accounts for less than 30% of the ther- 



mal pressure. But lSijacki et al.l (|2008[ ) affirmed that 
this component can reach up to 50% of the central 
gas pressure in clusters. 

Despite all the effort in computing the cos- 
mic ray pressure, the distribution of cosmic-rays 
in ICM is yet poorly known, and direct evidence 
for cosmic-ray ions in the ICM is still lacking. The 
measurements cited above indicate that cosmic ray 
pressure accounts fo r a minor contributi on to the 
dynamical support (jAndo fe Nagai 120081 ). We ex- 
pect that future experiments like the Imaging Air 
Cherenkov Telescopes (IACTfH which will work in 
the TeV band) and the Gamma-ray Large Area 
Space Telescope (glastE which will work in the 
GeV band) will be able to provide better constraints 
to P cr in clusters. 

Bearing these results in mind, we used central 
values for the ratio between the cosmic ray pres- 
sure and the thermal pressure between 10%-50%. 
To which end, we adopted 0.1 < Y p0 < 0.5. 

The value of depends on the model adopted 
for cosmic ray dynamics in clusters. In the sim- 
plest model \& = 0, the energy distribution of cos- 
mic rays follows precisely the thermal gas in the 
clu ster. The recent radiati ve simulations performed 
by IPfrommer et al. showed that ^ takes a 

value of -0.5. We considered \P = —0.5 based in 



( Pfrommer et al. 2007). 




2 http:/ /heasarc. gsfc.nasa.gov/docs/cgro/cossc/egret/ 

3 http:/ /magic. mppmu.mpg.de/introduction/iact. html 



http:/ /www-glast. stanford.edu/ 



'I'm 

Fig. 4. Variation in the mass profile due to the cos- 
mic ray pressure as the only non-thermal component 
for the clusters in our sample. The red line represents 
the hydrostatic mass profile, while the green line shows 
the maximum variation of the mass profile due to the 
contribution of the magnetic pressure. The blue zone 
represents the region of the mass profile variation de- 
pending on the ^ parameter. 



In Fig. [4] we show the mass profile variation due 
to the cosmic ray pressure as the only non-thermal 
component, compared to the hydrostatic mass pro- 
file. Comparing the results presented in Fig. [T] with 
Fig. 0] we verify that the cosmic ray pressure is the 
most important non-thermal component for NCC 
clusters (see also Tabled. 
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4. Mass determination including the effects 
of non-thermal pressure 

We have assumed that the cluster is spherically 
symmetric, the gas is in magnetohydrostatic equi- 
librium and that consequently, the gas pressure 
and the non-thermal components support the ICM 
against gravity. We can then write: 



d(P g +Pb+ Pturb + Per) 

dr 



-P9- 



GM NTP (r) 



(9) 



where the gas pressure (P g ) at a temperature T is 
PgkT/^rrip, the magnetic pressure (Pb) is {B 2 )/8tt, 
the turbulent pressure (Pturb) is \Pg{&1 +°f)' cos ~ 
mic ray pressure {Pgr) is Yp P g (r/r )*, G is the 
gravitational constant and Mntp is the total mass 
inside a radius r. In our case, we computed the to- 
tal masses within t^qq. Considering the influence 
of non-thermal pressures in the magnetohydrostatic 
equilibrium equation, we have the following expres- 
sion for the total mass of the cluster: 



Mp NT (r) 



k B T(r) f din p g {r) t dlnT(r) 



G/irriH 
dB(r) 2 



dlnr 

>,2 



8np g {r)G dr 2p g {r)G dr 



dlnr 



(p 9 (rK 2 (r)) 



G 



(2a 2 (r)-a t 2 (r)) 



r 2 dP cr (r) 
Gp g {r) dr 



(10) 



where p. is the mean molecular weight, m# is the 
hydrogen mass, fc is the Boltzmann constant, T(r) 
is the temperature profile, B(r) is the magnetic pro- 
file described in Sect. 13. ll and a r and a t are the ra- 
dial and tangential dispersion velocity of the intra- 
cluster gas, respectively (see Sect. I3.2[) . 

5. Results and discussion 

In Fig. [5] we show the evaluation of all non-thermal 
pressures separately in order to analyze each con- 
tribution alone (right panels). The maximum in- 
fluence which non-thermal components yield in the 
mass estimate is shown in left panels. From this 
Figure, we see that the main non-thermal contri- 
bution comes from magnetic fields or cosmic rays, 
depending on the range of parameters adopted. 

From this Figure we also see that in the in- 
ner parts (r <~0.5 r 500 ) of cool core clusters 
(A496, A1689 and A2667) the maximum cosmic 
ray pressure is higher than the maximum mag- 
netic pressure. On the other hand, for the outer 
parts of the mass profiles the maximum mag- 
netic pressure is always higher than the maxi- 
mum cosmic ray pressure. However, this statement 
is not true for non-cool core clusters. For A2050 



and A2631, the cosmic ray pressure dominates the 
magnetic pressure for all radii. Even the turbu- 
lent pressure becomes more important than the 
magnetic component (see Table [5]). It is impor- 
tant to note that we used results from the lit- 
erature to est imate the central magnetic field in- 
tensity 



<e- g - lTav 



or & Perlevi 119931: iFeretti et al. 
Il999at lAllen et all 120011: lEilek fc Owenl 12002^ and 
as observed by Vogt &: Enfilinl ( 20051 ) the strength 



of central magnetic fields in non-cool core clusters is 
lower than those present in cool core clusters. Thus, 
as the magnetic pressure is described by Eq. ([3|), a 
lower central value leads to a decrease in the mag- 
netic pressure. 

In Table we present the maximum differ- 
ence in mass estimates considering the influence of 
each non-thermal pressure separately and all three 
components together. We define erg (max) as the 
maximum difference in mass estimates due to the 
magnetic pressure only, <7t U rb(max) is the maxi- 
mum difference in mass estimates due to the turbu- 
lence pressure only, <7 cr (max) maximum difference 
in mass estimates due to the cosmic ray pressure 
only and Ototai (max) is the maximum difference in 
mass estimates due to all non-thermal components. 
The mass variation given by ct-Mntp(b) is simply 



a A/ntp = 



M NTP (r) -M(r) 
M(r) 



(11) 



Table 2. Maximum difference in mass estimates. 



Cluster 


ob (max) 


<7 turb (max) 


(j cr (max) 


OtotaJ ( max ) 


A496 


17.33% 


5.52% 


4.87% 


27.72% 


A1689 


20.07% 


3.87% 


9.47% 


33.40% 


A2050 


1.24% 


0.82% 


10.69% 


12.74% 


A2631 


0.59% 


0.77% 


9.79% 


11.15% 


A2667 


14.92% 


3.06% 


9.93% 


27.90% 



From Table [21 we see that for CC clusters the 
magnetic pressure is dominant, contributing more 
than 50% of the total mass variation. Still, for NCC 
clusters the major component is the cosmic ray 
pressure, accounting for more than 80% of the to- 
tal mass variation. We conclude that the maximum 
influence of non-thermal components on the total 
mass variation ranges from more than 10% to al- 
most 35%. Although all of the assumptions agree 
with previous works, we emphasize that the dif- 
ference in the hydrostatic mass estimates rely on 
the assumption of the specific parametrization of 
non-thermal components that were inferred based 
on numerical simulations and observational results. 
As the errors in mass estimates (by weak lensing 
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or via hydrostatic equilibrium) are in most cases 
lower than the variation in mass due to non-thermal 
pressure contribution, we cannot neglect this com- 
ponent. Moreover, as X-ray data are widely used to 
constrain cosmological parameters, this evaluation 
should be regarded with care. Since in most cases 
this method neglects the non-thermal contribution 
to mass estimates, these constraints may be biased 
low and consequently it will bias the cosmological 
determinations. Thus, multi- wavelength study will 
play an important hole in the investigation of non- 
thermal components in galaxy clusters. 

5.1. Constraining the non-thermal pressure 
based on A 1689 

The plasma in many apparently relaxed systems 
may be affected by additional non-equilibrium 
processes, which contribute to rise the pressure 
and hence cause an underestimate of the clus- 
ter mass from X-ray observations of the thermal 
bremsstrahlung emission. If this difference is due to 
non-thermal pressure, the comparison between in- 
dependent methods of mass estimates can provide a 
powerful constraint to the contribution of this com- 
ponent. 

It is important to notice that weak lensing study 
is only possible if the cluster is massive and not very 
close. For the sample chosen in this work, the only 
cluster with available weak-lensing mass determi- 
nation inside r 50 o is A1689. In this section we focus 
on the comparison of a hydrostatic X- ray mass esti- 



mate determined in a previous work (ILagana et al 



2008) with those derived from weak lensing which 
are available in the literature. However, as noted 
by iHoekstral (|2007l ). it is difficult to compare re- 
sults from different mass indicators since it would 
be necessary to make assumptions regarding the 
cluster geometry in the X-ray determination, and 
because lensing is sensitive to all mass along the 
line of sight, this would also bias weak lensing mass 
determination. 

The total mass inside rsoo computed from X-ray 
measurements is M 50Q = (11.14 ± 0.46) x 1O 14 M 
without co nsidering the non-th ermal component. 
For A1689. lMahdavi et all (|2008h estimated a total 
mass inside r 500 of M 500 = (14.29 ±2.40) x 10 14 M Q . 
In this case, the total non-thermal contribution 
(this means all three components: cosmic ray, mag- 
netic and turbulent pressure) to the mass estimates 
can range from ~2% to more than 30%. 

This cluster is possibly undergoing a merger, 
where a sub-clump close to or along the lin e- 



20041 ). Thus, in this specific case, the turbulent mo- 
tion seems to acc ount for the most p art of the non- 
thermal pressure. iNagai et alJ (|2007l ) affirmed that 
the hydrostatic estimate of the total mass is biased 
low by about 5%-20% through the virial region, pri- 
marily due to additional pressure support provided 
by subsonic bulk motions in the ICM. 

On average, the hydrostatic cluster mass esti- 
mates are biased low by about 7.5% at rsoo for 
relaxed systems, while the bias i n unrelaxed sys - 
tems is about 10.5% at this radius ( Lau et all 2009). 
From our results (see Fig. [5] for this cluster) we see 
that if one assume the non-thermal contribution for 
this cluster, the mass estimates from different meth- 
ods (e.g., weak lensing and X-rays) can be totally 
compatible. 

6. Conclusion 

We have taken into account the effects of non- 
thermal pressure on the X-ray mass estimates for 
five Abell clusters (A496, A2050, A1689, A2667 and 
A2631). The masses derived considering just the 
th ermal pressure were presented in a previous work 
bv lLagana et al. ( 2008 ) and were used here for com- 
parison. We summarize our main results below: 

— The inclusion of non-thermal pressure in the 
intra-cluster gas description is motivated by the 
increasing evidence for the presence of a mag- 
netic field in clusters of galaxies. We assume a 
magnetic profile given by B(r) cx B^p^ , consid- 
ering values for Bq ranging from 5 up to 30 /iG 
for CC clusters, while for NCC clusters we as- 
sume 2/iG < Bo < 8/iG. For each central value 
we let the shape parameter vary between 0.5 
< a < 0.9. The magnetic pressure contributes 
with approximately 20% to the total mass vari- 
ation. 

— In order to take into account the influence of tur- 
bulent motion in the ICM and bulk velocities, 
we assume isotropic turbulent press ure based 



on the numerical simulation results of lLau et al 
(2009): Pturb = l/3p g (cr, 2 + er t 2 )- Th e tangential 



of-sight is being accreted (jAndersson fc Madejski 



(er,.) and the radial (at) dispersion velocity pro- 
files were derived based in the same numerical 
simulations. This component can influence up 
to 5% of the total mass estimates. 
Energetic particles are confined by magnetic 
fields. Since ICM is permeated by magnetic 
fields, cosmic rays can also provide an impor- 
tant source of pressure. As the distribution of 
cosmic rays is poorly known, we cons idered the 
prescription of Ando &: Nagail (2008) to derive 
the cosmic ray pressure. From our results we 



Lagana, Souza & Keller: Non-thermal pressure in clusters of galaxies 



11 



can see that this component can affect the clus- 
ter mass estimates by ~ 10%. In the inner parts 
of cool core clusters this component is compa- 
rable to the magnetic pressure. 
The plasma in clusters of galaxies may be af- 
fected by non-thermal processes which rise the 
pressure and hence cause an underestimate in 
X-ray measurements of the total mass. In this 
way the comparison between independent meth- 
ods of mass estimates can constrain the contri- 
bution of non-thermal pressure. We compared 
weak-lensing results with X-ray mass measure- 
ments for A1689 to investigate this, and if the 
difference in mass estimates is in fact due to the 
non-thermal component, it can account for 2% 
to -30%. 

We took into account the effects of non-thermal 
pressure on the total mass estimates. To derive 
the thermal pressure and compute the hydro- 
static mass we used observed temperature and 



density profiles (jLagana et al.ll2008f ) . This is the 
first study that considers the influence of all 
three non-thermal pressure on the total mass 
estimates. Within the limits of our sample and 
with several reasonable assumptions to describe 
each of the non-thermal components, it is im- 
portant to notice that our findings rely on the 
specific parametrization adopted in this work. 
This study indicates that further investigations 
are needed to make a detailed description of the 
influence of non-thermal components in ICM. 
— Without the complete knowledge of the non- 
thermal contribution, we should be aware of the 
fact that clusters of galaxies have been used as 
observational tools for investigations of cosmo- 
logical interest based only on the hydrostatic as- 
sumption. The non-thermal components are ne- 
glected. A thorough knowledge of non-thermal 
contribution will come from the combination of 
X-ray, radio and gamma ray data for a large 
sample of clusters. The next generation of clus- 
ter surveys will provide data to address these 
fundamental questions in cosmologj{f|. 
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rAsoo r / r 500 

Fig. 5. Left panels: Minimum and maximum non-thermal pressure profiles for each component. The red dashed 
line represents the thermal pressure profile. For the non-thermal pressures, the blue continuous lines represent the 
maximum and the minimum cosmic ray pressure profiles, the thick black dashed lines represent the maximum 
and minimum magnetic pressure profiles and the dot-dashed green line represents the turbulent pressure profiles. 
Right panels: comparison between the hydrostatic mass profile (red lines) and the maximum mass profile due 
to non-thermal components. The blue region represents all mass profile variations due to all the combinations 
possible for the three non-thermal components. The clusters are displayed from top to bottom in the following 
sequence: A496, A2050, A1689, A2667 and A2631. 



